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We discuss dark energy models which might describe effectively the actual acceleration of the uni- 
verse. More precisely, for a 4-dimensional Friedmann-Lemaitre-Robertson- Walker (FLRW) universe 
we consider two situations: First of them, we model dark energy by phantom energy described by 
a perfect fluid satisfying the equation of state P = (/3 — l)p (with j3 < and constant). In this case 
the universe reaches a "Big Rip" independently of the spatial geometry of the FLRW universe. In 
the second situation, the dark energy is described by a phantom (generalized) Chaplygin gas which 
violates the dominant energy condition. Contrary to the previous case, for this material content a 
FLRW universe would never reach a "big rip" singularity (indeed, the geometry is asymptotically 
de Sitter). We also show how this dark energy model can be described in terms of scalar fields, 
corresponding to a minimally coupled scalar field, a Born-Infeld scalar field and a generalized Born- 
Infeld scalar field. Finally, we introduce a phenomenologically viable model where dark energy is 
described by a phantom generalized Chaplygin gas. 
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I. INTRODUCTION 

Several astronomical and cosmological observa- 
tions, ranging from the cosmic microwave background 
anisotropy to observation of distant supernova 0, 
show that the universe is undergoing an accelerating 
stage. In addition, these observations show that the 
acceleration of the universe is due to some unknown stuff 
usually dubbed dark energy (DE), which constitutes 
roughly two thirds of the total energy density of the 
universe. Moreover, it is known that the DE satisfies an 
equation of state P = ((3 — l)p, where \B\ < 0.3 (at least 
recently in the history of the universe) |3j . 

So far, several phenomenological models have been pro- 
posed to describe the dark energy, being the cosmological 
constant, A, by far the most simple and popular candi- 
date 0. However, this possibility is ruled out (in prin- 
ciple) as a consequence of the huge discrepancy between 
the expected theoretical and experimental value of A. 
A positive cosmological constant might describe the ac- 
celeration of the universe as it could be described as a 
perfect fluid with negative pressure, —A, and this is one 
of the main ingredients to produce an accelerating uni- 
verse. In fact, a Friedmann-Lemaitre-Robertson- Walker 
(FLRW) universe undergoes an accelerated stage as long 
as p + 3-P < 0, where p and P correspond, respectively, 
to the total energy density and pressure of the matter 
content. Matter contents with this requirement can be 
described effectively, for example, by a perfect fluid; with 



a barotropic equation of state or a more exotic one like 
in (generalized) Chaplygin gas models [U |(J , or by dy- 
namical scalar fields as in quintessence models Q and 
phantom energy models 1 [1 EGH E [H . 

If dark energy would be described by either of the last 
two models mentioned above, then the future of the uni- 
verse might be quite different. While for a quintessence 
scalar field, with an effective equation of state P = 
(/3 — l)p, with (3 constant and < (3 < 2/3, dominating 
the energy density of the universe, the universe would 
expand forever, for a phantom energy; i.e. (3 < 0, this 
might not be the case. In fact, for a matter content 
with a barotropic equation of state formally similar to 
the previous one, but with a negative (3, the universe 
would experiment a cosmic doomsday, also dubbed big 
rip, [I Efl [H E H El; i.e. the scale factor would 
blow up in a finite cosmic time. The last affirmation is 
based on a constant negative value of (3. However, the 
value of may change along the evolution of the universe 
and then in principle the universe might not reach a big 
rip in the future. 

Another candidate to describe dark energy is a (gener- 
alized) Chaplygin gas p|, already mentioned, which cor- 
responds to a perfect fluid with a rather strange equation 
of state P — —A/p a , where A is a positive constant and 
a a parameter. This fluid can describe a transition from 
a dust dominated universe at early time to a de Sitter 
universe at late time. In addition, this matter content 
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1 We would like to mention that there are oth er candidates for 
dark energy based on brane-world mod els Il5l and modified 4- 
dimensional Einstein-Hilbert actions ITtl . where a late time ac- 
celeration of the universe may be achieved. 
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has been proposed as a unification of dark matter and 
dark energy. In this paper, we will show that if the dark 
energy is modelled by a phantom generalized Chaplygin 
gas, then the universe will escape the big rip and will ex- 
pand forever. Others phantom energy models; i.e. matter 
contents with p + P < and a positive energy density, 
exhibiting similar property has already been proposed 
|l2t Il3j| . For example, this can be achieved considering 
a homogeneous minimally coupled scalar field with an 
appropriate potential Il2l or with a Born-Infeld homo- 
geneous scalar field However, in these models, the 
scalar field has the wrong kinetic energy. In this paper, 
we propose an alternative phenomcnological model to de- 
scribe phantom energy by means of a fluid which, firstly, 
satisfies the generalized Chaplygin gas equation of state 
and secondly, violates the dominant energy condi- 
tion |17| . Moreover, for this peculiar material content a 
FLRW universe would never reach a Big Rip. 

The paper can be outlined as follows. In the next sec- 
tion, on the one hand, we will review the phantom energy 
model with a constant equation of state {(3 is negative 
and constant) in FLRW universes, giving the explicit ex- 
pression of the scale factor for the three different spatial 
geometries. On the other hand, we will discuss if the 
presence of a positive cosmological constant might alle- 
viate the big rip problem. In section 11111 we introduce 
and study a dark energy model based on a generalized 
Chaplygin gas. In section IIVI we analyze our model in 
the light of scalar fields, corresponding to a minimally 
coupled scalar field, a Born-Infeld scalar field and a gen- 
eralized Born-Infeld scalar field. In section[V] we describe 
a phenomenologically viable model, where the dark en- 
ergy is described by a phantom generalized Chaplygin 
gas. Finally, we briefly summarize and discuss our re- 
sults in section IVll 



II. PHANTOM ENERGY 

Through the paper, we mainly consider the late time 
evolution of a homogeneous and isotropic universe. Fur- 
thermore, we model dark energy by phantom energy, 
which in the present section is described by a perfect fluid 
satisfying the equation of state P — ((3 — l)p, where (3 
is constant and negative. The conservation equation re- 
sults on p = Aa~ 3/3 , where A is an integration constant. 
Therefore, for /3 < the energy density grows with the 
scale factor instead of decreasing. For simplicity, we dis- 
regard the other matter contents of the universe as their 
energy densities decrease with the cosmic time and can 
be neglected in comparison with the energy density of 
phantom matter at very late time when a big rip could 
happen 2 . Consequently, the Friedmann equation can be 



expressed as 



2 In section |y] we consider the other matter components of the 
universe together with a phantom matter, defined in section 1*1111 
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where H is the Hubble constant and k = 1,-1,0, corre- 
sponding to spherical, hyperbolic or flat spatial sections 
of the FLRW model. 

For flat spatial section (k — 0) the scale factor scales 
with the cosmic time, t, as 



a(t) 



3/3/2 



to) 



2/(3/3) 



(2.2) 



where ao and to will be integration constants throughout 
the paper, corresponding to the initial radius and cosmic 
time of the universe 3 and C = (87rG/3)A As can be seen, 
for negative /?, the scale factor diverges in a finite cosmic 
time 

£oo0 = to — o/i^-n /o a o^ 2 ' (2-3) 



if t varies between to and iooo- Therefore, the universe 
would reach a cosmic doomsday [Toj| . For values of the 
cosmic time larger than iooo the scale factor will decrease 
until vanishing at t — > +oo. It can also be seen that when 
[3 approaches zero, t^o goes to infinity. A similar situa- 
tion can be found in the case of a FLRW universe with 
hyperbolic or spherical spatial geometry [see Fig.^. This 
is not surprising as for a vanishing (3, the FLRW geome- 
try corresponds to a de Sitter space-time sliced into flat 
sections and there is no longer a big rip. In addition, for 
a given initial value of the scale factor ao > exp[2/(3/3)], 
the larger is the value of /?, the larger is the value of i^o 
at which the big rip happens. Moreover, the latter the 
phantom energy starts dominating the energy density of 
the universe; i.e. the larger the value of ao, the sooner 
the universe reaches the big rip. 

In the spherical case (k = 1), the scale factor must 
be larger than a m i n = C 1 ^ 3 ^^ 2 ^. Otherwise, the Fried- 
mann equation (|2.1|l is not well defined. We have that 
the cosmic time scales with the scale factor as [20f 
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where F(6, c; d; e) is a hypergeometric series 20] . The 
cosmic time is finite whenever 4 — 1 < 1 — Ca 2_3/3 < 0, 



and constraints the model using the observational cosmological 
parameters. 

3 The integration constant to can be set equal to zero. However, 
ao must be different from zero, otherwise the scale factor will be 
vanishing at any cosmic time. 

4 A hypergeometric series F(6, c; d; e), also called a hypergeometric 
function, converges at any value e such that |e| < 1, whenever 
6 + c — d < 0. However, if < + c — d < 1 the series does not 
converge at e = 1. In addition, if 1 < o+c— d, the hypergeometric 
function blows up at |ej = 1 L'll . 



i.e. a m i n < a < (2/C) 1 ^ 2 ~ 3 ^ l \ For larger values of the 
scale factor the expression (|2.4|l breaks down and, conse- 
quently, we cannot immediately conclude either the ex- 
istence or the absence of a cosmic doomsday. However, 
the difference between the cosmic times corresponding, 
respectively, to ioo+ when the scale factor blows up, and 
to a given cosmic time t such that a(t) is larger than 
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can be expressed as follows 



zero and infinity in a finite cosmic time corresponding to 
1 1 1 



(71/(3/3-2) 



2 F n 



2 3/3-2' 3/3-2 
3/3 3/3 



3/3 V 2 ' 2(3/3- 2)' 2(3/3- 2) 



L-l 



1;-1 



(2.9) 



^00 + 



2(71/(3/3-2) 

3/3 
3/3 



3/3 



3/3-1 

3/3-2' 2(3/3-2)' 2(3/3-2) 



3|3 

" 2(3/3-2) 



1 



(2.5) 



In addition, it can be checked that the last expression 
is well defined, in particular the hypergeometric func- 
tion (see footnote^), whenever a is larger or equal than 

(2/C) 1 ^ 2 3 ' 3 ''. This value of the scale factor corresponds 
precisely to the maximum value allowed in Eq. I|2.4|l . 
Consequently, we can conclude that there is a big rip 
in a FLRW universe sliced into spherical sections filled 
with phantom matter when (3 is constant and negative. 
In this case, we have that the cosmic time elapsed since 
the scale factor acquires its minimum value a m i n at t = to 
up to the divergence of the radius of the universe is 
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Similarly, a universe filled with phantom matter with 
/3 constant reaches a big rip in the future, if the geometry 
corresponds to a FLRW universe sliced into hyperbolic 
sections (fc = — 1). In fact, on the one hand, we have 
that the cosmic time varies with the scale factor as 
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for a < C 1 ^ 313 2 \ On the other hand, we have that for 
larger value of the scale factor, the cosmic time reads 
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where ioo_ corresponds to the cosmic time when the scale 
factor reaches infinite values. The last two expressions 
are well defined at a = C 1 ^ 313-2 ^ (see footnote^}. In this 
model, it can be seen that the scale factor varies between 




FIG. 1: This figure shows the behavior of the cosmic time 
corresponding to the big rip as a function of the parameter 
j3 related to the ratio between the pressure and the energy 
density of the phantom matter. The dashed line corresponds 
to the case of a FLRW universe sliced into hyperbolic sections. 
The solid line corresponds to a homogeneous and isotropic 
space-time sliced into spherical sections. The cosmic time 
has been divided by C 1/{3f3 ~ 2) . 

Before concluding this section, we will analyze if the 
inclusion of a constant positive vacuum energy density; 
i.e. a positive cosmological constant A, may alleviate the 
big rip problem dues to phantom energy with constant 
equation of state (/3 constant). For simplicity and analyt- 
icity, we will restrict to the case of a FLRW universe with 
flat spatial geometry. The Friedmann equation reads 
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The solution to the last equation is 



a 3 * (t) = a^(l-DY 
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where A = A/3 and D is a positive constant given by 
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From Eq. I|2.11|l . it can be seen that the scale factor grows 
from an initial value ao and blows up in a finite cosmic 
time; i.e. the universe will face a cosmic doomsday, when 

t approaches i = t + (ln£))/(3/3VA)- For t < t, the 
scale factor decreases and the universe collapses when 
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t approaches infinite values. In addition, the larger is 
the value of ao, the smaller is t, and consequently, the 
sooner the cosmic doomsday happens. A similar conclu- 
sion holds for (3 (at least for ao > 1). Moreover, it can 
be checked that t approaches iooo, defined in Eq. I|2.3[) . 
when the cosmological constant vanishes. In summary, 
we have that the presence of a cosmological constant does 
not modify the general features of the model and the 
big rip cannot be avoided for j3 constant and negative. 
Moreover, a positive vacuum energy density cannot de- 
lay the happening of the big rip [see Fig. [2] . This can be 
explained as follows: the presence of a positive cosmo- 
logical constant in the model induces a bigger growth of 
the Hubble parameter and, consequently, the scale factor 
increases faster leading to a sooner big rip. 
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FIG. 2: This figure shows the behavior of the cosmic time 
corresponding to the big rip as a function of the cosmological 
constant. As the graphic shows, the largest value for t is 
achieved in the absence of a positive cosmological constant. 
Indeed, the larger is A, the sooner the big rip takes place. The 
cosmological constant and the cosmic time in the graphic are 
redefined as dimensionless quantities given by A/(Ca~ 3/3 ) and 

-3/3/2,7 



(3/2)/VCa P/ (t - to), respectively. 



III. GENERALIZED CHAPLYGIN GAS AND 
PHANTOM ENERGY 

The generalized Chaplygin gas can be described as a 
perfect fluid with the following equation of state 



P 



-A/ P a , 



(3.1) 



where A is a positive constant and a is a parameter. In 
the particular case a = 1, the equation of state (|3.1|) 
corresponds to a Chaplygin gas. The conservation of the 
energy momentum tensor implies 
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where ao and po are the initial scale factor and energy 
density, respectively. It can be checked that the dominant 
energy condition is fulfilled whenever A < p^ a+1 \ This 
requirement is strongly related to the initial values of 



the model and the specific equation of state through the 
constant 
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For positive values of B, P + p is positive and the dom- 
inant energy condition is satisfied. This is not the case, 
when B is negative. Let us see the behaviour of the en- 
ergy density for both cases. 

When the parameter B is positive, p will be a decreas- 
ing function of a. In fact, for 1 + a > the generalized 
Chaplygin gas interpolates between dust for small scale 
factors and a constant energy density at large scale fac- 
tors. This property has promoted the generalized Chap- 
lygin gas to be a candidate to unify dark energy and 
dark matter y|- For 1 + a < 0, the energy density be- 
haves on the other way round; i.e. p approaches A 1 /^^ 
for small scale factor and behaves as a pressureless fluid 
at late time. 

When the parameter B is negative, the energy density 
will be an increasing function of the scale factor. More- 
over, p is larger than A 1 ^ 1+a ' > when 1 + a < 0, reaching 
its minimum value at a = and blowing up when the 
scale factor approaches its maximum value 
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In what follows, we will disregard this set up (B < 
and 1 + a < 0). On the other hand, if 1 + a > and 
B < 0, the scale factor is larger than a, in such a way that 
p vanishes at this scale factor and approaches J 4 1 /( 1 + Q ) 
when the scale factor goes to infinity. We will henceforth 
analyze this last case, which can be included in the set 
of phantom energy models as p > and P + p < 0. 

As in the previous section, we consider a homogeneous 
and isotropic universe, where now the phantom energy 
is given by a generalized Chaplygin gas such that B < 
and — 1 < a. The Friedmann equation reads 
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If the FLRW universe is sliced into flat sections, then 
the cosmic time is related to the scale factor as 
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where D = 8ttG/3 and —1/2 < a. Firstly, we have that 
the scale factor is bigger than a defined in Eq. (|3.4|) . In 
this case there is no big rip: when the scale factor blows 
up, the cosmic time does too (see footnote |3J. In oppo- 
sition with the cases studied in the previous section, the 
present model does not show a cosmic doomsday because 
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the Hubble parameter approaches a constant non van- 
ishing value for large scale factors. Consequently, at late 
time the geometry of the model is asymptotically de Sit- 
ter. Although we have not been able to get an equivalent 
analytical expression to Eq. I|3.6|) for — 1 < a < —1/2, a 
similar conclusion holds because H 2 approaches a posi- 
tive non vanishing value when a — > +oo. 

When the spatial geometry of the homogeneous and 
isotropic space-time is spherical, the Hubble parameter 
is well defined as long as a > a m j n , where a m j n is such that 
H (a m i rl ) = 0. The explicit expression of a mm is given in 
the appendix. It can be shown that a m i n is larger than 
the minimum value of the scale factor for flat spatial 
geometry given in Eq. I|3.4jl . Moreover, the cosmic time 
for k — 1 satisfies the inequality 
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The second term on the right hand side (rhs) of the in- 
equality is finite. Indeed, it is the cosmic time for flat 
spatial sections corresponding to a = a m i n . In addition, 
the first term on the rhs corresponds to the cosmic time 
for k — geometry at a given scale factor a > a. As can 
be seen for large scale factor the cosmic time for k = 1 
blows up because t — 1 diverges for flat spatial geometry 
(first term on rhs). Consequently, the universe does not 
hit a cosmic doomsday in its future. 

Similarly, the cosmic time for a FLRW universe with 
spatial hyperbolic sections can be bounded from below 
as follows 
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when the matter content corresponds to a generalized 
Chaplygin gas. We would like to point out that the last 
two terms on rhs of the inequality coincides precisely with 
the ones on the rhs of the expression l|3.7|l . Consequently, 
based on an argument similar to that is used for the k = 1 
case, we can conclude that there is no big rip for k = — 1. 
In addition, the scale factor grows from a to infinity. 

Once analyzed the late time behaviour of a homo- 
geneous and isotropic space-time filled by a general- 
ized Chaplygin gas with the characteristics already men- 
tioned, let us see how it behaves for smaller scale factors. 



The energy density vanishes whenever the scale factor 
approaches a, which can be the case only for flat and 
hyperbolic sections. Consequently, at a = a the pressure 
may diverge inducing a singularity in the geometry [see 
Fig- EJ. Indeed this can be the case if a is positive. The 
scalar curvature for k — ±1,0 reads 
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where the dot represents derivative respect to the cosmic 
time. As can be seen, R is well defined at any scale factor 
for spherical geometry (we recall a, < a min ). On the other 
hand, for k = —1, 0, the scalar curvature R is finite (even 
at a = a) as long as — 1 < a < 0. The same can de 
deduced for positive values of a except at a — a, where 
there is a divergence of R. We would like also to point 
out that for flat sections the FLRW universe presents a 
bouncing at a — a, which is regular if — 1 < a < 0. 




FIG. 3: The behavior of the pressure of the generalized Chap- 
lygin gas is shown in terms of its energy density for negative 
values of the parameter B. The dashed line corresponds to 
a = —0.1. The solid darkest (lightest) line corresponds to 
a — 1 (a = 2). The pressure and energy density has been re- 
defined as dimensionless quantities given by PA- 1/(1+a) and 
pA~ 1 ^ 1+a \ respectively. As can be seen, for positive values 
of a the pressure reaches extremely negative value when p 
approaches zero. 

Before concluding this section, we analyze the parame- 
ter (3(a) = P/p+1, which somehow quantifies the devia- 
tion of the generalized Chaplygin gas from a cosmological 
constant [see Fig. 0] and can be expressed in terms of the 
scale factor as 



B 
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(3.10) 



As can be expected (3 is negative for the set of parameters 
we are considering. At late time, [3 approaches zero; i.e. 
the FLRW universe is asymptotically de Sitter. On the 
other hand, (3 blows up near a (only for k = 0,-1). This 
is partially due to our oversimplified model. In principle, 
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we should have considered the other matter contents of 
the universe, as dark matter component, which are the 
dominant components for smaller scale factors. Indeed, 
if we consider dark matter (DM) given by dust, the ef- 
fective 5 value of j3 approaches the unity; i.e. the total 
matter content behaves effectively as dust, when a — ► a, 
as long as — 1 < a < 0. For positive value of a the 
effective value of f3 is still divergent at a. This can be 
understood as a consequence of the divergence of P near 
a for positive a. 




1.1 1.2 1.3 1.4 1.5 1.6 1. 1 1.8 



FIG. 4: This plot shows the behavior of the parameter given 
in Eq. (13. 101 as a function of the dimensionless variable a /a. 
The lower graphic (the darkest one) corresponds to a — —0.1. 
The upper graphic (the lightest one) corresponds to a = 2. 
Finally, the middle graphic shows /3 for the Chaplygin gas; 
i.e. a — 1. As can be see, the generalized Chaplygin gas 
approaches a positive cosmological constant (f3 — 0) for the 
largest values of the scale factor. 



In summary, we have shown that a generalized Chap- 
lygin gas can describe phantom energy. In addition, this 
material content can avoid the occurrence of a cosmic 
doomsday in the future of the universe. This is not sur- 
prising as the energy density of a generalized Chaplygin 
gas approaches a constant positive value for B < and 
< a + 1 . Consequently, a FLRW universe filled with 
this gas is asymptotically de Sitter. 



IV. GENERALIZED CHAPLYGIN GAS AND 
SCALAR FIELDS 



Up to now, we have described the generalized Chap- 
lygin as a perfect fluid with a peculiar equation of state 
(|3.1|) . In the following, we will describe the generalized 
Chaplygin gas in terms of scalar fields. First, we will 
show how the generalized Chaplygin gas (with a negative 
parameter B, see Eq. (|3.3[l ) can emerge in the context of 
generalized Born-Infeld phantom theories. For this pur- 



5 The effective value of is denned as P/p + 1, where now p is 
the total energy density of the universe and P is the sum of the 
pressure of the different material components. 



pose, we consider the Lagrangian Ca defined as 
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where g^ v is the metric of the space-time and is a scalar 
field. For a FLRW universe, Ca reduces to 



1 



(<PY 



(4.2) 



The dot corresponds to derivative respect to the cosmic 
time. It can be shown that the ener gy d ensity ps and the 
pressure Pa associated to Ca reads [T|| , 
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Consequently, ps and Pa satisfy a generalized Chaplygin 
gas equation of state; i.e. Pa = —A/p^. The difference 
between the Lagrangian defined by expression 1)4. Q and 
the one given in [f| is that the kinetic energy term for 
the scalar field 4> is negative. In addition, it can be seen 
that 



PaIpa = -[i + {4>)- 



(4.4) 



This expression shows that the scalar field <f> behaves as 
phantom energy [see also Eq. (|4.5I) ]. Moreover, this char- 
acteristic of 4> allows negative values of the parameter B, 
defined in Eq. 1)3. 3|) . as has been considered in the last 
section. Additionally, on the one hand, the time deriva- 
tive of 6 scales with the scale factor as 
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On the other hand, for FLRW universes with flat or hy- 
perbolic spatial sections filled by a generalized Chaplygin 
gas with B < and — 1 < a, the scale factor can takes 
any value such that a < a (a is defined in Eq. (|3.4() ). 
Consequently, for — 1 < a < 0, <f> vanishes when a ap- 
proaches a. However, for positive values of a, the time 
derivative of <fi diverges when a approaches a. For large 
value of the scale factor the opposite behavior is found: 
i.e. (f> approaches zero when < a and diverges when 
— 1 < a < 0. The divergence of 4> is harmless for large 
values of the scale factor, as the geometry of the universe 
behaves like a de Sitter space-time and, consequently, 
there is no singularity. Additionally, in a FLRW uni- 
verse with a flat spatial geometry the scalar field (f> varies 
with the scale factor as 



A 2 (!+°) 



x F 



1 



1 



Vfl 



SV 3(l+ot) 



1 /a\ 3 ( 1+Q ) 
ITc? 2(1 + a)' 2(1 + a) Va 



(4.6) 
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where — 1 < a and 4>q is an integration constant corre- 
sponding to the value acquired by <p at a = a. The scalar 
field (j> is finite for any value of the scale factor, whenever 
a is positive (see footnote 0J. For —1 < a < 0, the last 
affirmation remains true except for very large values of 
the scale factor, where <f> blows up (see Fig.^J). 




FIG. 5: This figure shows the behavior of the scalar field 
4> as a function of the dimensionless scale factor a/a given 
in Eq. 1)4. 6jl . The graphic with dashed line corresponds to 
a = —0.1, while the graphic with darkest (lightest) full line 
corresponds to a = 1 (a = 2). In the plot, the scalar field p 
has been redefined as V^GA 1/Wl+a)) (0 - <f> ). 

In the following, we show how the generalized Chaply- 
gin gas can be described effectively in terms of a Born- 
Infeld phantom scalar field, ip, whose Lagrangian reads 



-VWy/l-g^V^Vvr/t. 



(4.7) 



For a homogenous and isotropic space-time, the energy 
density and the pressure associated to ip reads 



density 

mm 



v(4>) 
i + ip 2 



p 4 , = -v(ip)^i + ip 2 



(4.8) 



Obviously, for a Chaplygin gas; i.e. a = 1, with phan- 
tom energy characteristics, the potential V(ip) is con- 
stant; V — V~A (see Eq. Ij4.1|l for a = 1). However, in 
general a generalized Chaplygin gas can be described ef- 
fectively by a Born-Infeld phantom scalar field, ip, only 
when its potential depends explicitly on ip. In fact, it 
can be easily seen that V varies with the scale factor as 
follows 



V(a) =Va(A 



B 



a 3(l+a) 



(4.9) 



The potential V approaches a constant value for large 
values of the scale factor (we are considering — 1 < a and 
B < 0). In addition, its behaviour near the minimum 
value of the scale factor, a, (for k = 0,-1), depends 
strongly on the specific value of the parameter a: for 
\a\ < 1 the potential vanishes near a, for 1 < |a| the 
potential blows up (see Fig. 0. A main difference be- 
tween the behaviour of the scalar fields ip and ip is that ip 



always reaches large values near a, while this is not nec- 
essarily true for <j>. Additionally, it can be seen that for 
a FLRW universe with spatially flat sections, the scalar 
field ip varies with the scale factor as 



A 2 (! + o) 



\f6irGa 
a 



1 



^3(i+«) 



(;) 



2(1 + = ) 



F| I 1 • ° 

" 2' 2(1 + a)' 2(1 + a) 



8\3(l+«) 



(4.10) 



when a is positive. In the last equation ipo is a constant 
corresponding to the value acquired by at a = a. It can 
be shown that ip is finite for any value of the scale factor 
(see footnote^]). On the other hand, for — 1 < a < 0, the 
scalar field scales with a as 



Ip-lpo 



A 2(1 + Q ) 



VforG(l +a) 

2 + a 1 3 
2(1 + a)' 2' 2 



V a 



3(1+q) 



(4.11) 



where "0oo is the value reached by ip for very large scale 
factors. In this case, the scalar field ip is well behaved for 
any value of a, expect at a = a where it blows up. 




V[a] 



1.2 
1.15 

1 . 1 
1.05 
1 



12 3 



FIG. 6: These figures show the behavior of V(a) defined in 
Eq. 14.91 as a function of the scale factor. The figure on the 
top corresponds to the values a = —0.1 (full line) and a = 1/2 
(dashed line). The figure on the bottom corresponds to a — 2. 
The potential V(a) has been redefined as A" 1 ^ 1+a ^V{a) and 
the scale factor as a/ a. 

Finally, we analyze the behaviour of a phantom min- 
imally coupled scalar field, x, able to mimic the be- 
haviour of a generalized Chaplygin gas (for B negative 
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FIG. 7: These figures show the behavior of tp defined in 
Eqs. ipTTfljl and fTTTJl as a function of the scale factor for 
a FLRW universe with flat spatial sections. The figure on 
the top corresponds to a = —0.1. The figure on the bot- 
tom corresponds to a = 0.5 (dashed line) and a = 2 (full 
line). In the top plot, the scalar field tp has been redefined as 
\/6lrCL4 1/(2(1+Q,)) (-!/> - Vo). In the bottom plot tp is redefined 
as \/feG4 1/12(1+a,) (i|) - ^oo). In both plots the scale factor 
has been divided by a. 




and — 1 < a). In this case, the energy density p x and 
pressure P x of the homogeneous scalar field \ read 



Px = -\x 2 + V( X ), P x = -\x 2 -V{ X ). (4.12) 



If the scalar field x simulates a generalized Chaplygin 
gas, the potential V scales with the scale factor as 



V(a) = 



1 



A- 



B 



a 3(l+a) 



2A 



B 



a 3(l+a) 



.(4.13) 



As can be seen, V(a) approaches a constant value when a 
blows up. This is not surprising: as we have already men- 
tioned a FLRW universe filled with a generalized Chap- 
lygin gas is asymptotically de Sitter. For the scalar field, 
X, this results on V(a) approaching a non vanishing con- 
stant and a vanishing x for very large scale factors. In ad- 
dition, V(a) is finite when a approaches a (for k = 0, —1) 
as long as — 1 < a < 0. However, for positive values of a, 
the potential V blows up near a. Moreover, If the FLRW 
universe is spatially flat, \ is wen behaved for any value 
of a. In fact, its evolution can be described in term of 



FIG. 8: These figures show the behavior of V as a function of 
the scalar field if) [see Eqs. |JI3J, lETTtJl and lETTO ! for a FLRW 
universe with flat spatial sections. The figures starting from 
the top correspond to a = —0.1, a — 0.5 and a — 2, respec- 
tively. The potential V has been redefined as A^ 1 ^ 1+a ' l V. 
In the top plot tp is redefined as V / 6^GA 1/(2(1+Q,)) (tp - tp x ). 
On the other plots, the scalar field tp has been redefined as 
V^GA 1/{2 ^ +a)) (tp-tp ). 



the scale factor as 
X-Xo 



(1 + a)V2A^G 



arcsm 

2 



5 s 3(1+q) 



(!) 



1 



,(4.14) 



where xo is an integration constant. In a addition, the 
potential V varies with \ as 



V(X) 



1 fA\ 1 

2 \2 

3 — cos 



(1 + a)V2^G{ X - xo) 



{l - cos [(1 + a)V2A^G{x - Xo)] } ' 



(4.15) 
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The behaviour of the potential V(x) depends strongly on 
the value of the parameter a (see Fig-EJ • For — 1 < a < 




FIG. 9: These figures show the behaviour of V(x) defined 
in Eq. H4.15|l as a function of the minimally coupled scaler 
field x f° r a FLRW universe with flat spatial sections. The 
figure on the top corresponds to a — —0.1. The figure on the 
bottom corresponds to a = 2. In the plots, the scalar field 
X has been redefined as \/247rG(x — Xo) and the potential 
as 2(2/A) 1/<1+a) V(x). As can be seen the behaviour of V 
depends on the chosen value of a. For negative values of a 
the scalar field rolls up the potential. However, for positive 
values of a, the scalar field x rolls down the potential. 

the scalar field rolls up the potential. However, for posi- 
tive values of a, the scalar field \ rous down the potential 
in contrast with the result obtained in Ref. [IJ- In the 
first case, the scalar field starts with a vanishing veloc- 
ity (% = 0) climbing the potential. Its velocity continues 
increasing until it reaches a maximum value and then it 
starts decreasing, vanishing for very large values of the 
scale factor (see Fig. EJJ. In the second case, x starts 
with an infinite velocity rolling down the potential. Its 
velocity is continuously decreasing, until vanishing when 
the scale factor blows up (see Fig. II 01) . 



V. A PHENOMENOLOGICALLY VIABLE 
MODEL 

In order to study the possible occurrence of a big rip 
(in the future of the universe) caused by phantom energy, 
it is a good approximation to consider that the matter 
content of the universe is mainly given by phantom en- 
ergy at very late time (large scale factors). However, any 
cosmological viable model able to describe the actual ac- 
celeration of the universe has to take into account the 




1.5 



20 
15 



1. 1 



1.2 



1.3 



1.4 



1.5 



FIG. 10: The behaviour of the quadratic velocity of the scalar 
field ~x is shown as a function of a/ a. The figure on the top 
corresponds to a = —0.1, while the figure on the bottom 
corresponds to a — 2. x 2 has been redefined as A~ 1 ^ 1+a ^x 2 ■ 



other material components of the universe, in particular 
DE. Consequently, the Friedmann equation reads 



H z = 



8ttG 



(pde + Pdm + Pr) , 



(5.1) 



where pde,Pdm,Pr correspond, respectively, to the en- 
ergy density of DE, DM, and radiation. We will consider 
that DE is described by a generalized Chaplygin gas with 
B < and -1 < a (see Eq. (E2J)) and the DM compo- 
nent as a dust fluid. The Friedmann equation can be 
rewritten as 



H 2 = H? 



pc., \a/ \a 



(5.2) 



In the last expression, Hi, at, pc,i are the present Hub- 
ble parameter, scale factor and critical energy density. 
On the other hand, Odm and are the density pa- 
rameters for DM and radiation. The model can de- 
scribe the present acceleration of the universe as long 
as Pde/pc,4 — 0.7, flu. — 0.3 and £1r ~ 0. On the other 
hand, at the radiation dominated era (T = lMeV) the 
energy density must be dominated by the pr, (for a suc- 
cessful Nucleosynthesis). Considering that the present 
scale factor is equal to the unity, the scale factor in the 
radiation dominated era is equal to 2.4 x 10~ 10 . Our 
model can describe a viable cosmological model satisfy- 
ing all these requirements for different values of a [see 
Fig. , although there is a fine tuning of the parame- 
ters A and B. 
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FIG. 11: These figures shows the behavior of Ode, !2dm and 
Qr as a function of the scale factor since the nucleosynthesis 
epoch; i.e. p r = lMeV . The lightest graphics correspond 
to Q.r. The darkest full line graphics correspond to Odm • 
Finally, the dashed line graphics correspond to Qde- On the 
one hand, in the figure of the top, it has been chosen a — —0.1, 
~ 2.8 x IQ-^eV and |_B|V(i+«) ~ 3.8 x HT 40 eV. 
On the other hand, on the figure of the bottom, it has been 
chosen a = -0.9, A 1/{1+a) ~ 2.8 x lQ~ w eV and \B\ 1/{1+a) ~ 
1.4 x 10- 41 eV. As can be seen the behavior is very similar 
for both choices. 



The model involving the phantom generalized Chaply- 
gin gas can describe the actual acceleration of the uni- 
verse [see Sec.0, where the energy density of the gas cor- 
responds roughly to two thirds of the total energy density 
of the universe. However, there is a fine tuning between 
the parameters A and B related to the energy density of 
the gas. 

We have also shown how the phantom generalized 
Chaplygin gas can appear in a natural way in the con- 
text of phantom generalized Born-Infeld theories, where 
the kinetic energy term for the scalar field is negative. 
In addition, we have analyzed different effective phan- 
tom scalar fields which can mimic the behaviour of this 
gas. This has been carried out in the context of a Born- 
Infeld scalar field and a minimally coupled scalar field. 
Although, the dynamics of each of this scalar fields can 
be quite different, they provide the same behaviour for 
the scale factor and avoid the occurrence of a cosmic 
doomsday in the future of the universe. 

Finally, we would like to stress that a phantom en- 
ergy model does not necessarily imply a big rip in the 
future of the universe as has been shown using a phan- 
tom generalized Chaplygin gas, even without imposing 
any restriction on the sound speed |14| . Moreover, a big 
rip singularity is not necessarily related to a phantom 
matter content as has been recently pointed out in [19|. 
Although there the big rip singularity happens in a finite 
cosmic time, scale factor, energy density and Hubble con- 
stant; and the singularity is associated with a divergence 
of the pressure. 
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APPENDIX: EXPLICIT EXPRESSION OF a min 

The Friedmann equation (|3.5|) for k = 1 can be ex- 
pressed as H 2 = f(a) where 



f(a) = ~l/a 2 + D 



A 



B 



a 3(l+a) 



(A.l) 



The function /(a) has a unique positive root which we 
have denoted a m ; n . Its corresponds to the minimum ra- 
dius of a FLRW universe filled by a generalized Chaplygin 
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gas with B negative and — 1 < a. The explicit expression 
of <z m i n depends on the parameter j2]| 



Q = 



B 2 £)-3(l+«) 

lA? 27A 3 ' 



(A.2) 



For positive value of Q; i.e. (27 '/ '4)D 3 ( 1+Q ) AB 2 > 0, a 
reads 



B 
2A 



B 
2A 



For Q = that is £ 2 = 4/(27L> 3 ( 1+Q U), a min is 

1 / 4 \ 2(l + o) 



(A.3) 
(A.4) 



Finally, if (27 / '4) D^+^AB 2 < 0; i.e. negative value of 
Q, Omin can be expressed as 



2g 1/3 cos 



where 



V27£> 3 ( 1 +«)A 3 ' 



arctan a / — 1 



27 £)3(1+«) J 4 B 2 



(A.5) 



(A.6) 



It can be checked that for any value of Q, a m ; n is larger 
than a given in Eq. (|3.4|) . 
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